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ABSTRACT 

We reconstruct the evolution of r\ Car in the last two centuries, under the assumption that the 
two 19th century eruptions were triggered by periastron passages, and by that constrain the binary 
parameters. The beginning of the Lesser Eruption (LE) at the end of the 19th century occurred 
when the system was very close to periastron passage, suggesting that the secondary triggered the 
LE. We assume that the 1838-1858 Great Eruption (GE) was triggered by a periastron passage as 
well. We also assume that mass transferred from the primary to the secondary star accounts for the 
extra energy of the GE. With these assumptions we constrain the total mass of the binary system to 
be M = Mi + M2 > 250M Q . These higher than commonly used masses better match the observed 
luminosity with stellar evolutionary tracks. Including mass loss by the two stars and mass transfer 
from the primary to the secondary we obtain a good match of periastron passages to the two peaks in 
the light curve of the GE. Based on these finding and a similar behavior of P Cygni, we speculate that 
major LBV eruptions are triggered by stellar companions, and that in extreme cases a short duration 
event with a huge mass transfer rate can lead to a bright transient event on time scales of weeks to 
months (a "supernova impostor" ) . 

Subject headings: (stars:) binaries: general— stars: mass loss— stars: winds, outflows— stars: individual 
(77 Car) 



1. INTRODUCTION 

77 Car is a binary system (Damineli 1996) containing 
a very massive Luminous Blue Variable (LBV; hereafter 
'the primary') and a hotter and less luminous evolved 
main sequence companion (hereafter 'the secondary'). 
The LBV has undergone two major eruptions in the 19th 
century. The Great Eruption (GE) took place between 
1837.9 — ~1858 and created the bipolar Homunculus 
nebula which contains 10 — 40M Q , and possibly more. 
(Gomez et al. 2006, 2009; Smith et al. 2003; Smith 
& Ferland 2007; Smith & Owocki 2006). Smith (2009) 
has even suggested that the Homunculus can be more 
massive. Following the GE, the Lesser Eruption (LE) 
took place between 1887.3-1895.3. This was a much less 
energetic eruption (Humphreys et al. 1999) and only 
0.1 — 1M were ejected from the primary (Smith 2005). 
A summary of the observed visual magnitude from the 
earliest measurements in the 17th century, through the 
19th century eruptions and up to 2004 was prepared by 
Frew (2004). Presently r\ Car is being continually moni- 
tored in the optical wavebands by the La Plata observa- 
tory (Fernandez Lajus et al. 2010), and occasionally by 
other observatories. Figure [T] presents the visible obser- 
vations from the last two centuries. 

At present the orbital period is Pq = 2023 days = 
5.539 yrs (Damineli et al. 2008a; Fernandez Lajus et 
al. 2010; Landes & Fitzgerald 2009), and the orbit is 
very eccentric. We will use an eccentricity of eo = 0.9 
as in our previous papers (e.g., Kashi & Soker 2009). 
Counting back in time the periastron passages, simply by 
subtracting a natural multiplicity of Pq, it is evident that 
the LE started very close to periastron (Frew 2004). This 
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Fig. 1. — The historical V-mag light curve (Frew 2004, marked by 
black symbols; Humphreys et al. 1999 & Fernandez Lajus, private 
communication, marked by blue diamonds) during the last two 
centuries. The two major eruptions, the 1837.9 — ~1858 Great 
Eruption (GE), and the 1887.3-1895.3 Lesser Eruption (LE) are 
marked on the figure. The estimated errors of the observations 
of Frew (2004) are given in the legend. The two rapid rises in 
V-mag in 1837.9 and 1843 at the beginning of the GE are clearly 
seen. Even though that the error in the peaks is 0.4-0.5 mag, it is 
evident that the rise in these two peaks was > 1.5 mag. 

is shown by the dashed blue line in Fig. [5] The beginning 
of the LE just as periastron passage occurred, leads us 
to our basic assumption: The Major LBV eruptions of 
77 Car are triggered by the periastron passages of the 
secondary. The major eruptions are defined as eruptions 
in which the luminosity increases by a few magnitudes, 
as oppose to regular eruptions, e.g. S Dor phases, or 
weak eruptions on top of the GE light curve, in which 
the luminosity is changed by < 0.5 mag. The primary 
must already be in a very unstable phase in order for 
the secondary to trigger its major eruption. For that, 
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Fig. 2. — The binary separation (blue lines) and the V-mag light 
curve (Frew 2004, marked by black symbols, see legend in Fig. 
[T] Humphreys et al. 1999 & Fernandez Lajus, private commu- 
nication, marked by diamonds) during the Lesser Eruption (LE) 
of r) Car (1887.3-1895.3). Dashed blue line: counting periastron 
passages back in time, not taking into account any mass loss or 
mass transfer. We used the common parameters Pq = 5.539 yrs, 
o = 16.64 AU, e = 0.9, h'h = 120M Q and M 2 = 30M Q to draw 
the binary separation, but individually these parameters, apart 
from Pq, have no effect on the times of periastron passages. It is 
evident that the LE started very close to periastron (Frew 2004). 
Solid blue line: including mass loss (for a total of lMg) and mass 
transfer from the primary to the secondary (for a total of O.14M0) 
at a constant rate for the duration of the LE. The masses before the 
event are Mi = 170Mq and M2 = 80Mq, to match the values at 
the end of the GE in the model named 'MTz' described in section 
13.31 We assume that during the periastron passage the interaction 
of an already unstable primary star with the secondary triggered 
the LE. 



such a major eruption does not occur at each periastron 
passage. It is the companion that makes the eruption so 
energetic. 

We shall use the common definition that phase zero is 
dehned at periastron passage. Note that the so called 
spectroscopic event in present rj Car — the weakening or 
even disappearance of ionization lines and the minimum 
intensity in all bands — does not necessary occur exactly 
at periastron, but possibly a few days earlier or later 
(Damineli et al. 2008a, b). This has little importance for 
our model, in which we only interested when the perias- 
tron passages occur, and not in the spectroscopic event. 
We only assume that periastron took place a few days 
earlier or later than Jan 11, 2009, as observed (Fernan- 
dez Lajus et al. 2010; Landes & Fitzgerald 2009). As 
will be evident below, for our calculation the possible 
few days error in determining the exact time of present 
periastron passages presently, is much smaller than the 
other sources of error. 

Previously, under different considerations, Livio & 
Pringle (1998) speculated that the GE was triggered by 
strong interactions among the stars and their winds at 
periastron passage. Soker (2001, 2004, 2005, 2007) sug- 
gested that the GE was caused by disturbances in the 
outer boundary of the inner convective region, which ex- 
pelled the outer radiative zone. According to these pa- 
pers, this resulted in a mass loss of ~ 20M© from the 
primary star, of which 12M0 were accreted onto the sec- 
ondary that ejected 4M Q of them by blowing bipolar jets, 
and kept the other 8M Q . The accretion rate was high, 
and the potential lobe of the primary is probably filled 



for most of the orbit (the potential lobe is the Roche 
lobe analogue, as here there is no synchronization while 
the RLOF process assume synchronization). In addition 
Bondi-Hoyle accretion (wind accretion) was also present. 
We studied these two processes in the context of the 
present regular wind of rj Car, and found that they both 
exist close to periastron (Kashi & Soker 2009). The ac- 
cretion during the GE was probably an hybrid of Roche 
lobe overflow and Bondi-Hoyle accretion. The accreted 
mass released gravitational energy that might account for 
the entire extra energy of the GE, while the jets shaped 
the bipolar nebula around the star (the 'Homunculus') 
and supplied part of its kinetic energy. 

Damineli (1996) has already tried to connect periastron 
passages with the light curve of the GE, and marked the 
periastron passages on the light curve of the GE, accord- 
ing to a period of 5.52 yrs. He found that the periastron 
passages fall not far from rises in V-mag. However, he 
did not consider the mass loss process, and assumed a 
constant periodicity, and for that could not accurately 
reproduce the times of periastron passages. 

In this paper we go back in time to the GE and LE, 
and examine what were the binary parameters back then, 
under our basic assumption that binary companions trig- 
ger major LBV eruptions. We describe our calculation 
method and assumptions in section [2J In section [3] we 
present the results of our calculation. We discuss the 
implication of our results in sectional 

2. THE HISTORICAL ORBITAL PARAMETERS 

We start with the latest periastron passage in 2009.03 
(Corcoran 2009; Fernandez Lajus et al. 2010), and move 
back in time to the 19th century, taking a constant or- 
bital period of Pq — 2023 days = 5.539 yrs (Damineli 
et al. 2008a; Fernandez Lajus et al. 2010; Landes & 
Fitzgerald 2009). The effect of the LE on the binary pa- 
rameters is small, as only < IMq were ejected during 
this eruption (Smith 2005), and most likely much less 
then this amount was transferred from the primary to 
the secondary. Nevertheless we take the LE into consid- 
eration and check how it affected the orbital parameters. 
Going further back in time we reach the end of the GE, 
in 1858. As considerable mass was evidently lost from 
the system during the GE, and possibly transferred from 
the primary to the secondary (Soker 2004), the binary 
parameters have changed during the GE. As mass loss 
and mass transfer were continuous during the GE, so are 
the variations in the binary parameters. 

We perform a calculation to determine the change in 
the orbital parameters, as described below. We repeat it 
twice, one time for the LE and one time for the GE. The 
initial conditions for the calculation for the LE arc taken 
from present day system parameters. The calculations 
of the binary evolution during the GE are constrained 
by the demand that the binary parameters at the end 
of the GE are like those at the beginning of the LE as 
we derived first. Practically, we integrate the equations 
backward in time. 

Two major effects are responsible for changing the or- 
bital period and the other parameters (eccentricity and 
semimajor axis): mass loss from the binary system in- 
creases the orbital period, while mass transfer from the 
primary to the secondary reduces it. For the binary pa- 
rameters of rj Car, per unit mass, the effect of reducing 
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the orbital period by mass transfer is much stronger than 
the effect of increasing it by mass loss. 

To calculate the variation in the binary and orbital pa- 
rameters we follow the equations from Eggleton (2006). 
The rates of change of the stellar masses (going forward 
in time) are 



The specific angular momentum is 
h = r x r. 



Mi = -rhn - m t ; M 2 = -m l2 + m t 
M = Mi + M 2 = -mn - ma, 



(1) 



where rhu and mi 2 are the rates of mass loss to infinity 
from the primary and the secondary, respectively, and m t 
is the rate of mass transferred from the primary to the 
secondary. The mass lost from the secondary is assumed 
to come from the transferred mass. It is ejected by the 
accretion disk in a bipolar outflow. For clarity in our 
following explanations, we shall use the product of these 
parameters with the duration of the GE, Mu = mntcE 
and M t = to^ge- To minimize the parameters' freedom, 
we keep the mass loss rate of each star and the mass 
transfer rate constant during each of the 19th century 
eruptions. 

The orbital separation is calculated as a function of 
time. The orbital separation r varies according to 
(Eggleton 2006) 



GM(t)r(t) 



i 



i 



Mi{t) M 2 (t) 



r(t). (2) 



The first term on the right hand side is the unperturbed 
gravitational acceleration, and the second term is a per- 
turbing acceleration due to mass transfer and mass loss. 
The perturbing acceleration depends linearly on the ve- 
locity, and for that its effect is larger close to periastron 
passage. Because we know the orbital period and we have 
an estimate of the eccentricity at present, we performed 
the integration backward in time to just before the GE. 
In this procedure the present parameters of r\ Car serve 
as the initial conditions. As the changes in the masses 
are not negligible, the equation cannot be solved analyt- 
ically and it is solved numerically, using an automatic 
step-size Runge-Kutta-Fehlberg integration method, of 
order 2-3 (implemented by MATLAB function ODE23, 
using accuracy factor of 10~ 8 ). 
The eccentricity e(i) = |e(t)| is calculated according 

to 

GMr 

GM e = f x (r x f) . (3) 

r 

The Keplerian energy per unit reduced mass e{t) is cal- 
culated according to 



, , 1. 2 , N GM (t) 

and then we can calculate the semi-major axis 

GM(t) 



a(t) 



and the orbital period 



P(t) = 2tt, 



Mt) ' 



GM(t) ' 



(4) 



(5) 



(6) 



(7) 



We can summarize our method of calculation and as- 
sumptions as follows: (1) Mass transfer and mass loss 
rates are constant during the orbital period. This might 
not be exactly the case. Mass loss and mass transfer rates 
might have been larger near periastron passage, while the 
slow relative motion near apastron might have increased 
the mass accretion rate near apastron. (2) Mass trans- 
fer and mass loss rates are constant during the entire 
GE. This might not be exactly the case as it is likely 
that these rates decrease toward the end of the erup- 
tion. (3) The mass is rapidly lost. This is not exactly 
accurate near periastron passages, when the secondary 
orbital motion relative to the primary was larger than 
the equatorial wind speed. (4) The two very rapid rises 
by ~ 1.5 mag during the GE (see below) occurred dur- 
ing periastron passages. (5) Present day mass loss values 
are < 10" 3 M yr^ 1 (Hillier et al. 2001; Pittard & Cor- 
coran 2002), and therefore have negligible effect on the 
orbital parameters. In other words, present day param- 
eters are assumed to be valid at the end of the LE. (6) 
The extra energy of the two eruptions comes from mass 
accretion onto the secondary. This assumption helps to 
constrain the binary parameters, but the main result of 
this paper, that under the assumption that the peaks 
of the GE occurred near periastron passages we can fit 
the evolution of r/ Car, does not not depend on this as- 
sumption. This will lead us to conclude that very likely 
all major LBV eruptions are triggered by companions at 
periastron passages. The triggering of a major eruption 
requires the primary to be in an unstable state accord- 
ing to our model. That also results in very high mass 
transfer and mass loss rates. 

3. RESULTS 
3.1. The Lesser Eruption 

The mass transfer in the LE is ~ 0.04 times that in the 
GE, as the total extra energy was ~ 0.04 times that in the 
GE, if we use our assumption that the extra energy comes 
from accretion onto the secondary (see equation[§]below). 
Under that assumption the accreted mass during the LE 
is taken to be 0.14M Q . For the total mass ejected in the 
LE we take 1M Q (Smith 2005). The orbital separation 
as function of time during the LE, calculated with the 
prescription described in section [2j is plotted in Fig. [2j 
The parameters used are listed in the figure caption. The 
LE did not change the orbital period by much. For the 
parameters we used, we find that at the beginning of the 
LE (and at the end of the GE) the period was 5.52 yrs, 
compared with the present period of 5.54 yrs. 

When going back in time to the LE we find (follow- 
ing Frew 2004) that it began only ~ 1 month after the 
system was at the 1887.2 periastron. We assume this is 
not a coincidence but rather a strong hint that the peri- 
astron passage was the trigger of the LE. At periastron 
the secondary is very close to the primary, to the extent 
that the tidal force it exerts on the primary is able to 
reinforce mass transfer and accretion. We suggest that 
the periastron passage of 1887.2 occurred when the pri- 
mary was in a more unstable phase, e.g. strong magnetic 
activity (Harpaz & Soker 2009). The periastron passage 
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during this more unstable phase triggered the eruption. 

One word of caution is in place here. The light curve of 
the GE contains strong fluctuations in brightness (Frew 
2004). Namely, some of the peaks might be either obser- 
vational errors of stochastic variations in r\ Car. How- 
ever, the three peaks we based our analysis on (one in 
the LE and two in the GE) are much larger than the 
other peaks, and seem to be real properties of the erup- 
tions. Though for the two peaks of the GE the error is 
0.4-0.5 mag (Frew 2004), it is evident that the rise was 
> 1.5 mag. The error in the timing in the data by Frew 
(2004) is small to be of a significance for our results, 
which require accuracy of the order of ~ 0.1 yr. 

3.2. Models for the Great Eruption with no mass 
transfer and no accretion 

As the LE started close to periastron (Frew 2004), we 
are encouraged to think that a similar effect has triggered 
the GE. The GE has two peaks in V-mag with rise larger 
than 1.5 mag, at 1837.9 and 1843. We use our two ba- 
sic assumptions to constrain the masses of the two stars. 
The two basic assumptions are that (a) major LBV erup- 
tions are triggered by periastron passages, and (b) that 
most of the extra energy is released by mass transfer to 
the non-eruptive secondary. 

The GE changed the orbital parameters more dramati- 
cally than the LE. The reasons for that are (a) that more 
mass was lost to the Homunculus (Smith et al. 2003), 
and (b) the much brighter event resulted from higher 
mass transfer to the secondary according to one of our 
assumptions. We perform the calculation described in 
section [5] for different sets of parameters. First, we take 
the 'common model' to be the model with the conven- 
tional parameters Mi = 12OM and M 2 = 3OM (in 
our previous works we used these stellar masses, adopted 
from the literature), and neglect mass transfer. To be 
punctilious, the masses we take here are for present day 
rj Car, and are changed by less than 1M during the LE, 
so practically they can be referred to as the masses at the 
end of the GE as well for this model (which we rule out 
later). In the 'common model' a mass of M/i = 12M Q 
is ejected by the primary, and no mass is transferred 
to the secondary or lost from it: M12 = M t = 0. We 
find that the orbital period before the GE in this model 
was P1837.9 = 4.76 yrs, and that periastron occurred 
~ 8.5 months after the event had started. 

According to one of our assumptions, periastron pas- 
sages trigger large rapid rises in luminosity. For that the 
first orbital period during the GE should be ~ 5.1 yrs as 
the time period between the two rapid rises by > 1.5 mag. 
In addition, the peaks should be close to periastron pas- 
sages, i.e., the orbit should be in phase with luminosity 
peaks. If only mass loss from the primary is included, 
with no mass transfer, then we can fit the period but not 
the exact times of periastron passage. In that case the 
periastron passage from the model occurs after the lumi- 
nosity rise in 1837.9, which is not the desired outcome. 

Under our assumption that the mass loss rate is con- 
stant during the GE, the common model does not match 
observations. We could relax the assumption of a con- 
stant mass loss rate. However, more mass loss is expected 
during periastron passages, a process that enhances the 
increase in orbital period due to mass loss. This will re- 
quire an even shorter orbital period before the GE, with 
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Fig. 3. — The variation of the binary parameters (orbital pe- 
riod P, semi-major axis a, orbital separation r, eccentricity e, and 
specific angular momentum h) during the 20 year long Great Erup- 
tion of r\ Car, for models with no mass transfer from the primary 
to the secondary. Left panels: M = 15OM0 ('common model'), 
right panels: increasing the total mass to M = 250Mq . Under our 
assumption that a total mass of 12Mq lost in a constant rate for 
the duration of the GE, the periastron passages do not occur close 
enough to the two major rises in V-mag, in 1837.9 and 1843, as 
has occurred in the LE, for any of the cases with no mass trans- 
fer. Changing the assumption of a constant mass loss rate might 
improve the fit for the high mass model, but not for the common 
model. 

a degraded fit to the two luminosity peaks. 

Taking a larger mass in the Homunculus, and hence a 
higher mass loss during the GE, will also make the dis- 
crepancy between the common model and observations 
larger. We therefore examine a more massive system for 
the same mass loss and zero mass transfer. As evident 
from equation @, when rh t = the individual masses 
of the stars play no role, and the variation in the orbital 
parameters depends upon the total mass alone. We find 
that under the assumption of a constant mass loss rate 
and a total mass lost during the GE of 12M Q , a first 
orbital period of 5.1 yrs years during the GE is achieved 
for a pre-GE total mass of M = (Mi + M 2 ) ~ 250M o 
(right panels of Fig. [3]). But even for this case, the times 
of periastron passage do not match the times of the two 
rapid rises. As a criteria for a good fit we take the sum of 
differences between the times of periastron passages P\ 
and P2, and the times of the rises in V-mag, 1837.9 and 
1843.0, respectively. Our criteria for best fit is therefore 
a minimal value for the quality parameter 

AP = (Pi - 1837.9) + (P 2 - 1843.0) yrs. (8) 

We find that for the 'common model' AP = 1.12 yrs, and 
for a total system mass of 250M Q (with no mass trans- 
fer) AP = 0.92 yrs. These values are considerably large 
comparing to the cases with mass transfer, discussed be- 
low. We therefore conclude that in order to get a good 
fit for the GE mass transfer must be taken into account, 
and that the two stars of rj Car are much more massive 
than usually thought. 
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3.3. The basic mass accretion model 

We now turn to use our assumption that the extra 
energy of the GE, Ege, both radiation and kinetic en- 
ergy, is mainly due to accretion onto the secondary 
star (Soker 2001, 2004). Soker (2007) took a mass of 
30M© for the secondary star, and deduced that the 
accreted mass onto the secondary during the GE was 
M acc = M t — M12 ~ 8M , which is a significant frac- 
tion of the secondary mass. Our finding above that the 
binary system is much more massive than assumed be- 
fore, allows for a much lower total accreted mass onto 
the secondary, hence accretion caused small disturbance 
to the secondary. We expect, therefore, that the present 
secondary will have about its main sequence properties. 
Taking for the secondary luminosity L2 = 9 x 10 5 L Q 
(Hillier et al. 2001), we find its mass to be M 2 ~ 80M Q 
(e.g., Meynet & Maeder 2003). We can also take for the 
radius of the secondary, R2, its zero age main sequence 
(ZAMS) radius, and constrain the total mass accreted 
during the GE to 

M acc = M t - Mi 2 

- 3 7 f E ™ \ ( R * \ ( M * \ 1 M 
^8 x 10 49 erg/ Vl4.3R y V8OM y ' 

(9) 

Similar to Soker (2007), who assumed that ~ 2/3 of the 
material transferred is accreted onto the secondary, we 
assume that M acc ~ 0.65M t , and take M t = 5.7M Q and 
M12 = 2M Q . We note that for our basic model 'MTz', we 
use the most basic assumptions, one of them is that the 
secondary is a ZAMS star. This assumption is probably 
not accurate. In section 13.51 below we consider a more 
realistic model in which the secondary is an evolved main 
sequence star, rather than a ZAMS star. 

Let us comment on the high mass accretion rate onto 
the secondary, ~ 0.2 — O.3M yr -1 in our models. In re- 
cent years some other models are based on similar, and 
even higher, mass accretion rates onto main sequence 
stars. In particular we note the merger model for the 
eruption of V838 Mon. In the stellar merger model of 
Tylenda & Soker (2006) an 8M main sequence star ac- 
cretes ~ O.O5M in about 3 months. This is an accre- 
tion rate of ~ O.2M yr -1 , about equal to the one used 
here, but onto a star that is only ^0.1 times the mass 
of the secondary here. Therefore, it seems that the pro- 
cess where main sequence stars accrete at very high rates 
might occur, and should be studied in more details. 

We run a number of models, changing the value of 
Mi , to find which value best fits the two rises in V-mag, 
according to our criteria from equation . The fitting is 
shown in Fig. 2) Checking primary masses in the range 
Mi = 120 - 3OOM we find best fit for Mi ~ 17OM , 
with AP = 0.28 yr. By that we determine the value for 
Mi to complete our basic model 'MTz'. It is interesting 
to note that the difference between the stellar masses for 
this case is the same at Mi — M2 = 9OM as in the 
common values used. 

With the mass transferred to, and the mass accreted 
by, the secondary, with the masses of two stars, and with 
the orbital period at the end of the GE fixed, we only 
need to fix the total mass lost by the primary directly to 
the bipolar nebula, i.e., the Homunculus. We fit Mh = 
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Fig. 4. — Fitting the best value for Mi for fixed parameters 
Ah = 8OM0, M(i = 18M , M; 2 = 2M Q , and M t = 5.7M . 
AP is our criteria for a good fit, given in equation lpl)l. Checking 
primary masses in the range Mi = 120 — 300Mq we find best fit 
for Mi ~ 17OM , with AP = 0.28 yr. By that we determine the 
value for Mi for our basic model 'MTz'. 

2OM , well within the estimated range of the Homuncu- 
lus, and from that take the mass lost by the primary to 
the Homunculus to be M n = M H - M /2 = 18 M yr -1 . 
The specific angular momentum is decreasing in a con- 
stant rate during the eruptions. In the common model 
case where there is no mass loss the specific angular mo- 
mentum remains constant (the lower left panel in Fig. 
[5]); note that the value of the angular momentum h has 
a negligible change. For an interesting discussion of the 
perturbations to the orbital parameters in binary sys- 
tems see Matese & Whitmire (1984). For a discussion 
including the occurrence of non-conservative mass trans- 
fer, and its consequences on parameters of binary sys- 
tems, see Sepinsky et al. (2009, 2010). 

We calculate the orbital parameters of the binary sys- 
tem using these parameters back to the beginning of the 
GE using the prescription described in section [2] The 
model with the above parameters is termed 'mass trans- 
fer ZAMS model' (MTz), and it is presented in Fig. [5] 
and[71 and summarized in Table[T]together with the other 
models. 

We emphasize that 'MTz' is not our preferred model. 
It simply served as the basic model to which other mod- 
els will be compared to. The model 'MTz' will serve to 
check the role of the different parameters below, and it 
will serve as the basis to examine more realistic mod- 
els that include the expectation that the two stars al- 
ready evolved from the ZAMS. In this paper our pre- 
ferred model is 'MTe' (see below), and we intend to ex- 
amine more models in future papers. 

It is clearly seen from Fig. [7] that the two sharp rises 
by ^ 1-5 mag in luminosity occurred very close to peri- 
astron passage. There is no perfect match, but it should 
be recalled that we did not perform too much of a pa- 
rameter fitting. For example, we tried to fit masses in 
increments of only 1OM , and only tried two values for 
the ejected mass. We only make sure that the scaling 
of the parameters as described above lead to an almost 
perfect fit. We could make further adjustments to the 
mass transfer and mass loss by the secondary, as well as 
with the stellar masses. Small variations in the param- 
eters and small changes from a constant mass loss and 
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A COMPARISON BETWEEN THE MODELS USED IN THE PAPER. 



mass transfer rates can lead to a perfect match. 

The main finding here is that a very massive binary 
system, where the primary experienced heavy mass loss 
and mass transfer, and where the secondary lost part 
of the transferred mass, fits the observations very nicely 
under our assumptions. The accreted mass onto the sec- 
ondary can account for the extra energy released during 
the GE. Further more, mass transfer and accretion onto 
the secondary make the fitting of periastron passages to 
observed luminosity peaks much easier than do models 
that include no mass transfer and accretion. 

3.4. The role of the different parameters 

We turn to analyzing the behavior of the orbital pa- 
rameters, orbital period P, eccentricity e, and semi- 
major axis a, during the GE. Their variation for our basic 
model 'MTz' are presented in the middle column of Fig. 
[5l together with the orbital separation r. The variation 
of the orbital period P during both eruptions is presented 
in Fig. |U We also study the role of the binary masses 
and amount of transferred and lost masses. For our ba- 
sic parameters ('MTz' model), a unit of mass transferred 
from the primary to the secondary reduces the orbital 
period much more than a unit of mass lost from the sys- 
tem increases it. Namely, M t plays a greater role than 
Mn and M12 in determining the orbital period. 

As the parameters space is large, we restrict ourselves 
to comparing the quality of the fit by changing one pa- 
rameter at a time. We therefore checked four different 
sets of parameters, to illustrate their influence on the 
results. The four additional models we checked are: 

1. Model 'Comparison 1' ('CI'): Decreasing the pri- 
mary mass to Mi = 120M Q . This results in a 
shorter binary period before the GE, as the mass 
transfer effect of decreasing the orbital period is 
smaller (proportional to M-f 1 — M% ; more nega- 
tive value implies a larger effect), and the mass loss 
effect that increases orbital period is larger. 

2. Model 'C2': Decreasing the secondary mass to 
M2 = 3OM . This results in a longer binary pe- 
riod before the GE, as the mass transfer effect of 
decreasing the orbital period is larger. 

3. Model 'C3': Decreasing the mass lost by the pri- 
mary to Mn = 12M Q . This results in a longer 
binary period before the GE, as in this model mass 
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Fig. 5. — The variation of the binary parameters (orbital period 
P, semi-major axis a, orbital separation r, eccentricity e, and spe- 
cific angular momentum h) during the 20 year long Great Eruption 
of r) Car. The variations are given for different sets of parameters 
we use in the paper (see Table[TJ. Left: 'Common model', middle: 
'MTz' model, right: 'MTe' model. 

transfer that acts to decrease the orbital period is 
more important. 

4. Model 'C4': Decreasing the mass transferred from 
the primary to the secondary to Af acc = 0M Q 
(namely, no mass transfer, only mass loss). This 
results in a shorter binary period before the GE, as 
no mass transfer implies that the mass loss effect 
is more important, and the latter works to increase 
the orbital period. 

Table Q] gives the parameters for the four cases, together 
with the value of the quality of fit parameter, AP (equa- 
tion[5]). As the parameters space is large, but it is evident 
that the parameters set for 'MTz' gives better fit then 
the other sets of parameters, our solution is probably a 
local minimum for the quality of fit parameter. Within 
the possible ranges of parameters, it seems as this is the 
solution for the system and mass loss and transfer pa- 
rameters. 

Other parameters, fixed by observations, are the be- 
ginning and termination time of the GE. As the orbit is 
very eccentric, the number of periastron passages play 
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Fig. 6. — The historical change of the orbital period for the 'mass 
transfer ZAMS model' (MTz). The present day orbital period Po = 
2023 days = 5.539 yrs is valid, when going back in time up to the 
end of the LE (1895.3). During the LE it is slightly reduced, but 
the more considerable reduction, still going back in time, occurred 
during the GE. The "~" like features occur close to periastron 
passages, and are more evident during the GE, when the mass 
transfer rate is higher than the mass transfer rate during the LE 
(both assumed to be constants) (see text for full explanation). 

a major role. The GE started very close to periastron, 
and lasted for more than 3 revolutions, but less than 4. 
As according to our model it started shortly after peri- 
astron, this means that there were 3 periastron passages 
during the GE. The most steep change in the orbital pe- 
riod occurred close to periastron passages. This is not 
surprising, since, as we mention above, the perturbing 
force (per unit mass) in equation [5J linearly depends on 
the orbital velocity, which is considerably larger close to 
periastron in such an eccentric orbit. 

We note that the value of the erj itself does not have any 
effect on the average change of the orbital period for a 
whole number of revolutions. The high eccentricity only 
causes the period to change at a different rate over the pe- 
riod. Taking any value of eccentricity does not affect the 
times of periastron passages. We only need an eccentric 
orbit for the secondary to be much closer to the primary 
during periastron, so it can trigger the eruption, but any 
value between eo = 0.8 and eo = 0.95 is acceptable for 
our model. The upper limit comes from the requirement 
that the orbital separation at periastron would be suf- 
ficiently smaller than the primary radius. The primary 
radius is ~ 0.83 AU today, taking primary luminosity of 
Li ~ 4.5 x 10 6 L Q (Davidson & Humphreys 1997) and 
effective temperature of ~ 20 000 K, and was probably 
larger during the GE. The lower limit comes from the 
requirement that the secondary would be much closer to 
the primary during periastron, to result in a significant 
tidal effect to trigger the eruption. Note that the same 
range of eccentricity comes from fitting doppler shifts for 
lines in present day rj Car (Kashi & Soker 2008b and 
references therein). The relative influence of the mass 
loss and mass transfer varies with orbital phase. For 
both processes the change of the orbital parameters per 
unit time (assuming constant mass loss and mass trans- 
fer rates) near periastron is larger than near apastron. 
The difference between periastron and apastron is larger 
for the mass transfer process. Therefore, in the first cy- 
cles of the GE the mass loss dominates near apastron 
in increasing the orbital period, while the mass trans- 
fer dominates near periastron and the orbital period de- 



creases for a short time near periastron (upper middle 
panel of Fig. [5]). This results in the "~" like features in 
the plot of the orbital period vs. time (Fig. [6]). As evi- 
dent from the middle column of Fig. [5j the magnitudes 
of the orbital period and semimajor axis decrease near 
periastron diminish over the cycles. This is explained as 
follows. The effect of the mass loss process on increasing 
the orbital period is proportional to (Mi + Ma) -1 , while 
that of the mass transfer in decreasing it is proportional 
to M-f 1 — M2 . As the total mass decreases during the 
GE the effect of the mass loss increases, while the effect 
of mass transfer decreases because Mi decreases and M^ 
increases. This causes the "~" like features which occur 
close to periastron passages to weaken from periastron 
to periastron, up to disappearing in the last one. 

The different rates in the changes (sometimes of oppo- 
site sign) of the orbital parameters between periastron 
and apastron make an exact fitting difficult. The reason 
is that fitting is sensitive to whether the GE ended near 
periastron or near apastron. Namely, the fit is sensitive 
to the way the mass loss rate and mass transfer rate de- 
clined at the end of the GE. For example, a whole num- 
ber of periods (for constant mass transfer and loss rates) 
would lead to no change in the eccentricity. However, 
the number of periods during the GE was ^3.5, and the 
eccentricity changed, but by a tiny amount for the 'MTz' 
model. It varies by the maximum value of Se = 5 x 10~ 4 
relative to its present value e = 0.9. However, having 
a large eccentricity the semi-major axis acquired most of 
its change in jumps close to periastron passages. 

If the GE had terminated earlier (later) , the final or- 
bital period would have been shorter (longer). The ec- 
centricity, however is less sensitive to the exact termina- 
tion time of the GE for the parameters we used, as its 
amplitude of change is small compared to its value. As 
we divide the mass loss and transfer in the duration of 
the GE, the duration only determine their rate, but not 
how much in total mass was lost or transferred. Accord- 
ing to our model 'MTz' the secondary was in apastron 
at 1856.88, and the following periastron was in 1859.63. 
The exact end date of the GE is not very important as 
long as it does not include another periastron passage, 
namely, as long it is before ~ 1859.5. As the end date of 
the GE is commonly considered to be as late as 1856-8, 
we conclude that the exact end date has a minor effect 
on our results. 

We conclude that it is not possible, under our assump- 
tions, to use a low mass system of ~ 150M Q to fit the 
rises in the light curve of the GE. Also, for our basic 
model 'MTz' with a total binary mass of 25OM , we find 
that mass transfer from the primary to the secondary is 
a necessary process to fit these rises. 

3.5. Our preferred model: an evolved secondary star 

The model 'MTz' serves as the basis model, for com- 
parison with more realistic models. In model 'MTz' (sec- 
tion I3.3[) the secondary star is taken to be a ZAMS star 
with a radius of 14.3R0. We now consider the evolved 
nature of rj Car, and consider an evolved secondary star, 
as we know the primary is already in its LBV stage of 
evolution and must be evolved. Based on the results of 
Verner et al. (2005) we take the secondary radius to be 
23.6R , but keep its post-GE mass at M 2 = 80M Q . To 
account for the energy released in the accretion onto the 




Fig. 7. — The binary separation (right axis; for two models) and 
the V-mag light curve (left axis; data from Frew 2004, see legend 
in Fig. [TJ during the Great Eruption of r) Car for the 'mass trans- 
fer ZAMS model' (MTz) (solid blue line), and the 'MTe' model 
(dashed gray line). The parameters are given in the legend in the 
order (Mi, M 2 , M n , M l2 , M t ); see Table Q] for definitions of sym- 
bols. 'MTe' is our preferred model, taking into account that the 
two stars already evolved from the ZAMS. For our preferred model, 
as well as the 'MTz model', two periastron passages nicely fit the 
rapid rises in V-mag in 1837.9 and 1843. 

secondary, the accreted mass should be M acc = 6.1M Q 
([9]). This value is more in line with the values suggested 
by Soker (2007). 

As mass transfer acts to increase the orbital period, a 
higher accreted and transferred masses should be com- 
pensated for by a change in other parameters, such that 
an increase in the primary mass and in the total mass lost 
to the Homunculus. We make two additional changes in 
the evolved 'MTe' model, from the basic 'MTz' model. 
(1) We take more mass to be lost to the Homunculus, 
as implied by the results of Smith & Ferland (2007) and 
Gomez et al. (2009). We emphasize that we do not know 
the Homunculus mass, but rather fit it (but our calcula- 
tions favor a large Homunculus mass). (2) We increase 
the primary mass. The increase in the primary mass 
is required in light of the very massive Homunculus in 
this model. As we see, the pre-GE mass of the primary 
is ~ 250M Q . Therefore, in our preferred model, 77 Car 
turns out to be a very massive binary system. 

As before, we do not look for a perfect match, and 
change Mi by increments of 10M Q . We find a good 
fit of the orbital passages to the GE peaks, as shown 
in Fig. for the following parameters: A primary 
and a secondary post-GE masses of Mi — 200M Q and 
M2 = 80M Q , respectively, where the primary has lost 
Mn = 38M© directly to the Homunculus and transferred 
M t = 8.1M Q to the secondary, out of which the sec- 
ondary accreted M acc — 6.1M Q (equation [5]), and lost 
the rest M l2 = 2M Q . 

The results of this section strengthen those of the pre- 
vious section. Models with an evolved massive binary 
system, a massive Homunculus, and where the extra en- 
ergy of the GE is supplied in a large part by accretion 
onto the secondary star, can account for the two peaks 
of the GE. 

4. SUMMARY AND DISCUSSION 



The Lesser Eruption (LE) of ij Car started close to 
the periastron passage of 1887.3 (Frew 2004), as evident 
from Fig. [2j This finding suggests that the interaction 
with the secondary, the less massive star, near periastron 
passage triggered the major eruption. We generalize this 
conclusion, and assume that the two rises in luminos- 
ity of the Great Eruption (GE) in 1837.9 and 1843 were 
also triggered by periastron passages. In addition to the 
strong tidal interaction exerted on the primary near peri- 
astron passage, major eruptions such as the LE and GE 
require the primary to be in a very unstable phase. 

In addition to the assumption that the GE peaks in 
the light curve occurred near periastron passages, we as- 
sume for simplicity that mass was lost from the system 
and that mass was transferred to the secondary at con- 
stant rates during the GE (to make parameters fitting 
minimal), that the total mass in the bipolar nebula (the 
Homunculus) is 12-40M Q (Smith 2005; Smith & Ferland 
2007; Gomez et al. 2009), and that the mass transferred 
onto the secondary released gravitational energy that can 
account for most of the extra energy of the GE (Soker 
2007). Not only the energy budget, but also the shape 
of the light curve, supports a mass transfer source for 
the extra energy. We note, though, that the assump- 
tion of accretion is not a necessary one for the results 
of the present paper, and it mainly used to constrain 
the amount of mass transferred from the primary to the 
secondary during the GE. 

The 1837.9 peak increased sharply and ~ 4 months 
later declined. A single star model is unlikely to pro- 
duce this very narrow first peak, as typical pulsation of 
unstable star do not have this type of behavior. Accre- 
tion onto the secondary, on the other hand, can easily 
account for that because the duration and shape can be 
accounted for by an accretion event (Soker 2007). After 
the first peak, the high and continuous mass loss rate 
from the LBV initiated close to the next periastron pas- 
sage of 1843. The following two periastron passages may 
have prolonged the eruption, possibly by extending the 
period of high mass loss rate. 

Solving for the orbital period evolution as a result of 
mass loss by the two stars and mass transfer from the pri- 
mary to the secondary, which effectively gives a shorter 
orbital period before the GE (section [3]), we could fit the 
1837.9 and 1843 sharp rises in the V-mag light curve with 
the occurrence of periastron passages (Fig. [7]). However, 
in order to achieve that fit we had to use stellar masses 
much larger than the previously used values. 

We used two models. Model 'MTz' assumes (wrongly) 
that the secondary is not evolved, and we took it to be 
a ZAMS star. This is our basic model; other models 
here are compared to it, as well as models we intend to 
run in future papers. In our preferred and more realistic 
model 'MTe', the secondary is evolved, with a radius and 
temperature as in Verner et al. (2005). We find that a 
very good fit for the motel 'MTz' is obtained for post- 
GE stellar masses of Mi ~ 170M©, M 2 ~ 80M Q , mass 
of 2OM lost in the GE, and 5.7M Q transferred from the 
primary to the secondary. For model 'MTe' a good fit 
is obtained for Mi ~ 200M Q , Mi ~ 80M Q , with 40M Q 
lost in the GE, and 8.1Mq transferred from the primary 
to the secondary. A detailed comparison between the 
parameters of the models is given in Table [TJ 

In the calculations we neglected tidal friction as a re- 
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suit of tides the two stars rise on each other close to 
periastron, as well as the drag force from the surround- 
ing ejecta that is not accreted. Both effects act to reduce 
the orbital period. However, as we estimate below, their 
effect is not expected to be large, although it can be 
non-negligible in the sense that the required transferred 
mass might be reduced by ~ 5 — 30% if these effects are 
included. The drag force due to particles that are not ac- 
creted is given by Alexander et al. (1976), and we used 
it as it was used by Livio & Soker (1984) 

F d = «ln(^)R% n vl indl (l--^Y, (10) 

V -^BH / \ ^2,Edd / 

where Rbh is the Bondi-Hoyle accretion radius, i? max 
is the cut-off distance, i> w indi is the relative velocity be- 
tween the secondary and the primary wind, L 2iC is the 
secondary luminosity, and £2,Edd is the secondary Ed- 
dington luminosity. It is assumed that gas particles 
within a radius Rbk are accreted by the secondary. The 
drag force has the role of Macc^windi that is played by 
the accreted mass in reducing the orbital separation. 

The large uncertainty is in the cut-off distance -R max . 
This is the distance beyond which the gravity of the sec- 
ondary is not the main force that determines the trajec- 
tory of the gas. Traditionally, -R max is taken to be the 
mean radius of the secondary's potential lobe. However, 
it seems that during the GE i? max was small, and the 
drag force weak. The cause for that was the huge lumi- 
nosity of the system, whether it comes from the primary 
or the secondary. 

More specifically, in our model the secondary has 
a super-Eddington luminosity, explained by accretion 
through the equatorial plane and an accretion disk. Most 
of the radiation escapes in a wide angle over the sec- 
ondary's polar directions. Indeed, most of the mass of the 
Homunculus is in two polar caps (Smith 2006). There- 
fore, the mass that was lost by the primary and passed 
near the secondary, was mainly flowing in the equatorial 
plane, and was accreted. Most of the influence on the 
orbit occurs near periastron, namely, at small radii. The 
high accretion rate implies that the effective accretion 
radius (-Rbh) was very large. Had the material been left 
around the secondary, it would have been able to influ- 
ence its orbit by exerting tidal forces. However, as the 
accretion radius is large the material was not left around 
the secondary but rather accreted onto it, and therefore 
it was not left there to exert the tidal forces. 

During the GE the binary system was within an 
optically thick region, and the luminosity was super- 
Eddington. This implies that radiation pressure on the 
gas in the polar directions was significant. The radiation 
and jets that were launched by the accreting secondary 
were therefore effective in pushing the gas in the polar 
directions, further reducing the amount of gas that is in- 
fluenced by the gravity of the secondary. Over all, we 
estimate that i? max was not much larger than -Rbh, and 
that i n (JW) < i. 

We also note that it is impossible to replace accretion 
by drag alone. The gravity of the secondary must influ- 
ence the flow up to a large distance, in order to make the 
drag force significant without accretion taking place. In 
that case we would expect accretion to occur. It seems 



one cannot escape the conclusion that significant accre- 
tion must have occurred during the GE. 

For reasonable parameters under our assumptions, we 
constrained the stellar masses to post-GE values of M\ ~ 
150 - 2OOM , and M 2 ~ 70 - 90M Q . We note that 
a secondary with a ZAMS mass of M% ~ 70 — 85M 
will have the required luminosity we use here of L 2 — 
9 x 1O 5 L (Verner et al. 2005). This is also the upper 
luminosity limit in the regime plotted for the secondary 
in the H-R diagram (Mehner et al. 2010), though this 
authors found a lower favorable value of L 2 — 4 x 1O 5 L . 
For a primary luminosity of L\ ~ 4.5 x 1O 6 L (Davidson 
& Humphreys 1997), a ZAMS mass M Y ~ 230M Q is 
required (Figer et al. 1998). We note that Figer et al. 
(1998) considered for the Pistol star a mass of ~ 200 — 
25OM and a metallicity slightly above solar metallicity. 
It seems as if very massive stars with solar metallicity 
can be formed in our Galaxy. Figer (2005) posed an 
upper limit for the most massive single star at 15OM . 
One possibility for forming such massive primary is a 
coalescence of two very massive stars, in a collision or 
a merger process. A similar process was discussed by 
Tutukov & Fedorova (2008) as a possible mechanism for 
forming massive blue stars. The high rotational velocity 
of the primary, ~ 0.8 — 0.9 of critical velocity (van Boekel, 
et al. 2003; Aerts et al. 2004; Smith et al. 2004) is 
also supported by this scenario. It is interesting to note 
that recently Crowther et al. (2010) observed a few very 
massive stars in the R136 star cluster, with the most 
massive member having presently a mass of ~ 265M , 
but possibly started as a 32OM ZMAS star. 

The largest uncertainty is the variation of mass loss 
and mass transfer during the orbit. The effects of both 
mass loss and mass transfer are much larger if occur 
near periastron passages. As we have not much informa- 
tion, we simply used a constant values during the GE. If, 
for example, tidal interaction near periastron passages 
caused burst of mass ejection by the primary, but the 
accretion onto the secondary continues for the entire or- 
bit, then we can do with less total mass loss and less 
massive primary. A resembling process occurs in cata- 
clysmic variables, where mass transfer rate depends on 
the shape of the Roche lobe of the donor, and sensitive to 
many other parameter, resulting in a varying mass loss 
rate (Ritter 1988). Because of the large uncertainties 
in the variations along the orbit of the mass loss and 
mass the transfer rates, we presently cannot strongly 
constrain the stellar masses. We can only say with a 
high confidence that before the GE the total mass was 
M x + M 2 > 25OM . 

With the parameters from the GE for model 'MTz' 
at hand, we reexamined the fitting of the LE. We could 
match the beginning of the LE with a periastron passage 
(solid blue line in Fig. [2]). The same fitting was per- 
formed for model 'MTe'. This shows that the LE is not 
sensitive to the uncertainties of the mass transfer and 
mass loss during the GE, and that our assumption that 
periastron passages trigger major eruptions is a reason- 
able assumption. 

The reason we find a high secondary mass comes di- 
rectly from our assumption that the luminosity of the GE 
was a result of mass accretion. For that, the larger the 
mass of the secondary the lower the required accreted 
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mass and vice versa (equation |HJ) . Mass transfer works 
to decrease the orbital period, while mass loss works to 
increase it. If we take the mass transfer to be too large it 
will not allow a good fit to the binary period we assume 
existed at the beginning of the GE (with periastron pas- 
sages close to 1837.9 and 1843.0). For that our model 
requires that the mass of the secondary be larger than 
the value usually assumed, in order to account for the 
GE luminosity. 

A very massive binary systems is supported by, or com- 
patible with, other considerations: 

1. Traditionally the masses of the two stars are de- 
termined by the assumption that each star shines 
at its Eddington luminosity limit (Davidson & 
Humphreys 1997). This assumption, however, is 
not in agreement when considering the evolution- 
ary paths of main sequence stars. Instead, we took 
the approach that the stars possess their luminos- 
ity as given by stellar evolution schemes. For a 
secondary stellar luminosity of L 2 — 9 x 10 5 L© 
(Verner et al. 2005), the MS secondary mass should 
be M 2 ~ 85M (Meynet & Maeder 2003). The 
stellar evolution calculations for very massive stars 
conducted by Yungelson et al. (2008; see also 
Figer et al. 1998) show that a star with a main 
sequence mass of somewhat more than 200M Q 
would become an LBV with the same luminosity 
of L\ ~ 4.5 x 10 6 L and effective temperature 
T\ ~ 20,000 K as that thought for the primary of 
■q Car (Verner et al. 2005). 

2. A more massive primary would make it easier to 
account for the super-Eddington luminosity dur- 
ing the GE, Lqe — 2 x 10 7 L©, because the ratio 
LoE/iEdd is smaller. 

3. In Kashi & Soker (2009) we showed that the mas- 
sive binary system model (termed there 'high- 
masses model'), favors accretion onto the sec- 
ondary close to periastron passage in present time. 
Such an accretion can explain the termination of 
the secondary wind that is required (Soker 2005a; 
Kashi & Soker 2009) to account for the deep de- 
cline in the X-ray emission (Corcoran et al. 2001; 
Corcoran 2005). 

We speculate that major LBV eruptions (also termed 
'giant eruptions'; Smith & Owocki 2006), are the result 
of unstable LBVs that are perturbed by the interaction 



with a companion during a periastron passage. This does 
not happen every periastron passage because the primary 
LBV must enter an unstable phase for that to occur. 
Our suggestion requires further study before it can be 
confirmed. Very few detailed major LBV eruptions have 
been observed. More observations of LBV binaries are a 
necessity to study this type of eruptions. For example, 
the series of major eruptions of the LBV in NGC 3432 
(Pastorello et al. 2010), which is possibly also a binary 
system, may be used to learn more about these eruptions. 

Because of the relatively detailed documentation of the 
GE, it serves as a prototype for major LBV eruption. We 
suggest that the LE of ij Car and the 17th century erup- 
tion of P Cygni (Kashi 2010) were also triggered by a 
periastron passage of the companion during a very un- 
stable phase of the primary star. In extreme cases the 
companion can accrete a considerable amount of mass 
and liberate a huge amount of gravitational energy dur- 
ing a short time, forming an optical transient event. In 
Kashi et al. (2010) we suggest that a somewhat similar 
process induced the optical transient NGC 300 2008OT- 
1 (observed by Monard 2008; Bond et al. 2009), and can 
account for other optical transients as well. 

According to our calculation, and under our assump- 
tions X] Car is the most massive binary system in the 
galaxy, leaving behind other massive stars, such as the 
Pistol star (Figer et al. 2004; Najarro 2005), WR 102ka 
(Barniske et al. 2008) and LBV 1806-20 (Figer et al. 
1998). However, we also predict that all other mas- 
sive stars that went through a non-spherical high mass 
loss rate phase have massive companions, e.g., P Cygni 
(Kashi 2010). Nota et al. (1995) already suggested that 
the non-spherical nebulae around LBV stars are shaped 
by binary interaction. So, it might still be possible for 
the Pistol star, for example, to take the lead. 
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